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Abstract 

It is shown, that the saturated a-effect taken from the non- 
linear dynamo equations for the thin disk can still produce 
exponentially growing magnetic field in the case, when 
this field does not feed back on the a. For negative dy- 
namo number (stationary regime) stability is defined by 
the structure of the spectra of the linear problem for the 
positive dynamo numbers. Stability condition for the os- 
cillatory solution (positive dynamo number) is also ob- 
tained and related to the phase shift of the original mag- 
netic field, which produced saturated a and magnetic field 
in the kinematic regime. Results can be used for expla- 
nation of the similar effect observed in the shell models 
simulations as well in the 3D dynamo models in the plane 
layer and sphere. 



1 Introduction 

It is believed, that variety of the magnetic fields observed 
in astrophysics and techni cs can be explained in terms o f 
the dynamo theory, e.g. ( HoUerbach & RiidigerJ 2004 ). 
The main idea is that kinetic energy of the conductive mo- 
tions is transformed into the energy of the magnetic field. 
Magnetic field generation is the threshold phenomenon: 
it starts when magnetic Reynolds number Rm reaches its 
critical value R^. After that magnetic field grows expo- 
nentially up to the moment, when it already can feed back 
on the flow. This influence does not come to the simple 
suppression of the motions and reducing of R^, rather to 
the change of the spectra of the fields closely connected to 
constraints ca used by conservation of the magnetic energy 
and helicity (B randenburg & Subramanian.ll2005h . The 



other important point is effects of the phase shift and co- 
herence of the physical fields before and after onset of 
quenching discussed in ( Tilgner & Brandenburg J20"08 ). 



As a result, even after quenching the saturated velocity 
field is still large enough, so that R^ » R^. Moreover, 
velocity field taken from the nonlinear problem (when 
the exponential growth of the magnetic field stopped) can 
still generate exponentially growing magnetic field pro- 
viding that feed back of the magnetic field on the flow is 
omitt e d (kinema t ic dyn a mo regim e) (Catt aneo & Tobias] 



20091, 'Tilaner,' bOOSt iTilgner & Brandenburg.! l2008t 



Schrinner, Schmi dt CameronJl2009l) . In other words, the 
problem of stability of the full dynamo equations includ- 
ing induction equation, the Navier-Stokes equation with 
the Lorentz force differs from the stability problem of the 
single induction equation with the given saturated veloc- 
ity field taken from the full dynamo solution: stability of 
the first problem does not provide stability of the second 
one. 

Here we consider effect of such kind of stability on an 
example of the model of galactic dynamo in the thin disk, 
as well as some applications to the dynamo in the sphere. 



2 Dynamo in the thin disk 

One of the simplest galactic dynamo models 
is a one-dimensional model in the thin disk 
dRuzmaikin. Shukurov. Sokolofi^.ll988l) : 



dA 

— =aB+ A 
dt 



dB 
dt 



(1) 



— = -DA' + B" 



where A and B are azimuthal components of the vector 
potential and magnetic field, a{z) is a kinetic helicity, D 
is a dynamo number, which is a product of the amplitudes 
of the a- and w-efFects and primes denote derivatives with 
respect to a cylindrical polar coordinate z- Equation ([T]i 
is solved in the interval -1 < z < 1 with the boundary 
conditions B - Q and A' = at z = +1. We look for a 
solution of the form 



(A, B) = e-^'iMz), S{z)). 



(2) 



Substituting (|2]i in ([T]i yields the following eigenvalue 
problem: 

(3) 

yS = -2)^' + S", 



1 



where the constant y is the growth rate. So as a(-z) = 
-a{z) is odd function of z, the generation equations have 
an important property: s ystem Q is in variant under trans- 
formation z — > -z when (|Parkerll97lh : 



j?[(-z) = j?[(z), S(-z) = 



or 



(4) 



J?[(-z) = -^(z), S(-z) = S(z). 



Therefore, all solutions may be divided into two groups: 
odd on S(z), dipole (D), and even, quadrupole on S(z). 
Then we can replace -1 < z < 1 with the interval < 
z < 1 and the following boundary conditions at z = 0: 
^' = 0, S = (D) and J?l = 0, S' = (Q ). Usually, a ^ 



ao with ao(z) = sin(7rz) is used, see also (lSoward.ll 19781) 
for ffoCz) = z dependence, more appropriate for analytical 
applications. 

System (O has growing solution, ^Ry > 0, when |2)| > 
For 2) < the first exciting mode is quadrupole 
with Tf' ~ -8 and 3y - 0: solution is non-oscillatorjQ. 
For 2) > the leading mode is oscillatory dipole, 
with higher threshold of generation: Tf^ ~ 200. Putting 
nonlinearity of the form 



a(z) = 



1 + 



for \B\ » 1 



(5) 



in ([U, where £■„, = (Z?^ + A'^)/2 is a mag- 
netic energy, gives stationary solutions for Q-kind 
of symmetry and quasi-stationary solutions for 
D, see about various forms of nonlinea rit ies i n 
dBeck. Brandenburg. Moss. Shukurov. Sokoloff. Il996l) . 



The property of the nonlinear solution is mostly defined 
by the form of the first eigenfunction. 

Now, in the spirit of ( Cattaneo & TobiasJ 2009 



Tilgner & Brandenburgj2008 ) we add to ([T]) equations for 
the new magnetic field (A, E) with the same a (|5]), which 



depends on (A, B) and does not depend on (A, B): 

dA 



dt 



^aB+ A" 



dB 

dt 



— = -OA' + B" 



M 
dt 



(6) 



— =aB+ A'' 



dB 

Yt 



— = -DA' + B" 



Numerical simulations demonstrate that for the negative 
D the both (A, B) and (A, B) are steady, however the 
final magnitudes of (A, B) depend on the initial condi- 
tions for (A, B), see Fig. [T] The procedure was the fol- 
lowing: equations for (A, B) were integrated up to 
the moment t - to, then the full system (|6]l was simu- 
lated with the initial conditions for (A, B) in the form: 
(A, B)\ = (A, B)\ (1 + Cb), where e e [-0.5, 0.5] 

l'='o l'='o 

is a random variable and C is a constant. The both vec- 
tors (A, B) and (A, B) are stable in time, however the final 
magnitude of E,„ for C ^ slightly depends on C. Pres- 
ence of alignment of the fields (A, B) and (A, B) follows 
from linearity and homogeneity of equations for (A, B), 
where a(z, E„) is given. Latter we consider stability of 
(A, B) in more details. 

For D > situation is dif jFerent, resembling 
that o n e of inst a bility d escribed in ( Cattaneo & Tobias.l 
2009 1: 'Tilgner,' '2008"; 'Tilgner & Brandenburg.l 120081: 



For our Galaxy usual estimate is D -- 



-10. 



Schr inner. Schmidt. Cameron. 2009) for more sophisti- 
cated models: field (A, B) oscillates and starts to grow 
exponentially, see Fig.|2] Note, that no regime in oscilla- 
tions for (A, B) is observed. The other specific feature is 

delay of (A, B) relative to (A, B): ^ 

If Em in Q is averaged over the space, so that a is 
steady, then instability dissapears. The question arises: 
does instability depends on stationarity, either it depends 
on something else? 

It is known, that for 3D < stability of the system 
(|3] |5]l, which has stationary solution, is tightly bound 
to behaviour of the linear solution of O) for £) > 
(iReshetnvak. Sokoloff. Shukurov.1ll992b . Note, that for 
the complex form of (HJ it is equivalent to the solution 
of the conjugate problem. 
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Figure 1: Evolution of magnetic energy E,„ for t < to 
governed by the system ( I1I5I) for £) = -10. In the moment 
fo - 300 the new magnetic field (A, B) with the initial 
conditions defined by constant C is switched on, see ( I6l5l l. 
Plots for E„, and E„ (with C = 0) for f > fo coincide. All 
the solutions are stationary. 



Let (A, B) = + fl, S + b), where S) is a solu- 
tion of the nonlinear problem and (a, b) is a perturbation 
with the same boundary conditions as for S). Putting 

— — Q(y 

(A, fi) in (O with o- ^ q-q + — yields equations for 
ja — a^b + a" 

(7) 

yb = -Dfl' + b". 



where a - a + ——S for a = 

dS 1 + S 



(1 + S2) 



:«() 



— T for 

S2 



|S| » 1. 



(8) 



Behaviour of o-w-dynamo (O is defined by the sign of Da, 
and its change in the perturbed equations ^ is important. 
In other words, instead of nonlinear equations (I3l5l l we 
come to the Unear problem (|3]l with given a = a(z, E„) 

D 

and effective dynamo number If - Then stabil- 

ity of fields (A, B) for the negative D can be explained 
as follows. For negative D solution (A, B) is finite and 
stable, because the threshold of generation for ^ is 
much larger than If, D^' » £)^ Field (A, B) is defined 



Figure 2: Evolution of magnetic energy £,„ for f < fo gov- 
erned by the system (1115b for D - 300. In the moment 
fo - 30 the new magnetic field (A, B) with the initial con- 
ditions defined by constant C is switched on, see (I6l5l l. 
Plots for £■„, and (with C = 0) for f > fo coincide. 
For C after intermediate regime the phase shift be- 
tween Em and E,,, increased and exponential growth of E„ 
started. 



up to an arbitrary factor, what corresponds to alignment of 
the vectors (A, B) and (A, B). Note, that D"^ «: does 
not guarantee, that (A, B) will grow exponentially due to 
nonlinearity (|5]l. 

It is worthy of note that nonlinear solution of ( 1115b and 
(1615b demonstrates similar stationary behaviour even for 

D 10^ in spite of the fact, that for the quadrupole 

oscillatory mode for positive D is ~ 200. The reason is 
that dynamo system tends to the state of the strong mag- 
leaving U 



netic filed with S ~ D^'^, so that a 



at 



waves. Here, instead of (|8]l we get 



J . If phase 



the level of the first mode's threshold of generation. 

For positive T) (A, B), and therefore a(B'), oscillate and 
one needs additional information on correlation of the 

shift between B and B is negligible, then a-effect is sat- 
urated and time evolution of (A, B) and (A, B) is similar. 
However, simulations demonstrate Fig.|2] that field (A, B) 
delays relative to (A, E). This is typical situation, when 
parameter resonance takes place: a is modulated by signal 
with frequency Q ~ 2<y, uj - 3y, see dDawes & Proctor 



-It is usually supposed that in aw-dynamo models B » A'. 



l2008h for details of spatial resonance. This assumption is 
supported by the fact, that instability disappears when in 
(|5]l steady a, averaged in time, is used. Note, that usage of 



3 



quadrupole boundary conditions in ^ is not important for 
instability: problem ( 1615b with periodical boundary con- 
ditions has oscillatory solution and instability depends on 
the form of quenching in the same way. 

To demonstrate what happens we consider how delay 6 
of (A, B) relative to (A, B) changes production of A^ + 
near the threshold of generation 2D". We start from the 
linear analysis of the system in the form: 



field (A, B). 



For D 



< (A, B) is steady and (A, B) is 
' » 1. 



unstable when , 

When (A, B) oscillates then (A, B) continues to oscil- 
late with (A, B) increasing the phase shift between (A, B) 
and (A, B). Then instability caused by the parameter res- 
onanse may arise. 



ia)A 
icjjB 



aB - k^A 
-ilTkA- k^B. 



(9) 



From condition of solvability for (|9]): (k^ + ioSf - 

—Hrf^kaQ with a - cq follows that ap- - k'^ - \. The 

other prediction of the linear analysis is the phase shift ip 

-~ — n 
between A and B: w - ±— , what is twice smaller than for 

i-^l — II 1 

the nonlinear regime jTilgner & BrandenburgJI2008l) . so 



that for the nonlinear regime the maximal A is when B is 
zero and quenching is absent. 

Then putting in ^ B - b sin(x - t), A - sin(x - t + 

1 



(fi + 0), B - sin(x - t + 9), and a = 



we get how 



generation depends on 6. Equation for B does not include 
original field (A, B), so we consider only production of 



2n 



A^. Then SAiip, 9) ^ [ dt. If |n| » 1, where 

_ 1 +B 

M(cp, 9) 

n = J , then (A, B) is unstable. 
SA(ip, 0) _ 

The exact equation for 6 A is: 

M((p, 9) - hi + /i2tan(^), 

cos(0)2(4 - 321/2) _ 2(21/2 _ 



h 



ho = 



21/2(21/2-1) 

sin(26')(32'/2-4) 

23/2(21/2 _ 1) • 



(10) 



If 61 = then hi = and A{(p, 0) = hi = 21^2 _ 4 jhen 



, h2 = -^(2''^ - 1), n at 



for = - hi = 
3 

ip +—is singular and instability appears. 

Summarizing results for the steady and oscillatory dy- 
namos we have the following predictions for stability of 



3 Conclusions 

Here we argue, that stability of the kinematic a<y-dynamo 
problem with the or-effect taken from the the weakly- 
nonlinear regime near the threshold of generation can be 
predicted from the knowledge on the threshold of gener- 
ation of the linear problem with the opposite sign of the 
dynamo number. It appears, that in spite of the fact, that 
the magnetic field akeady saturated a, it still can gener- 
ate magnetic field if spectra of linear problem are similar 
for dynamo number D with the opposite sign. So, as D 
depends on the product of the a and <y eff'ects similar anal- 
ysis can be performed with the c u-quen ching, usually used 
in geodynamo models, see e.g. dSoward, 1978.) . as well as 
with the feed back of the magnetic field on difiiision. It is 
likely, that for the more complex systems, velocity field, 
taken from the saturated regime, with many exited modes 
will always generate magnetic field if the Lorentz force 
would be omitted. 

So as nonlinearity ^ has quite a general form, we con- 
sider applications of these results to some other dynamo 
models. 

Linear analysis of the axi-symmetric al aw-equations 
gives the following, see (lMoff'att.lll978l) and references 
therein: for positive D (which is believed to be in the 
Earth) in presence of the meridional velocity Up the first 
exciting mode is dipole with 3y - 0. Reducing of Up 
firstly leads to oscillat ory dipole so lution (regime of the 
frequent reversals (.Braginskri.ll 1 964l) ) . The further reduce 
of Up gives the quadrupole oscillatory regime with larger 
value of 2If\ For negative D and Up the first mode 
is quadrupole with 3y = 0. t/p — » gives non-oscillatory 
dipole mode with decreased D", see for more d etails 



(iMeunier. Proctor. Sokolofi". Soward. Tobias.l Il997h . In 
contrast to the dynamo in the disk the thresholds of gen- 
eration for positive and negative D in the sphere are of 
the same order and situation with stability of the field B 



4 



is uncertain, and can depend on the particular form of the 
a- and w-effects. Anyway, stability of B for the steady 

regime is more expected^ 

In accordance with dCattaneo & Tobias.l l2009b shell 



models of turbulence demonstrate exponential growth of 
the magnetic field. This case, as well as 3D simulations of 
the turbulence in the box, which have the same instabili- 
ties, correspond to the oscillatory regimes and using our 
predictions should be unstable. 

In the case of the 3D dynamo in the sphere simula- 
tions demonstrate different behaviou r of B ( Tilgner.12008 



Schrinner. Schmidt. Cameron ] l2009h . For small Rayleigh 
numbers, when the preferred solutions is dipole (in oscil- 
lati ons) and close to the single mode structure (Case 1 



in ([Christensen. Aubert. Cardin. Dormv. Gibbons et al.. 
200 ihl B is finite. Increase of t he Rossby number 



(ISchrinner. Schmidt. Cameron.ll2009h leads to the turbu- 
lent state and B becomes unstable. 

Author is grateful to A.Brandenburg and D.Sokoloff for 
discussions. 
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